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The photogalvanic effect is studied in electron gas over the liquid He surface with the presence 
of quantizing magnetic field. The gas is affected by the weak alternating microwave electric field 
tilted towards the surface normal. Both linear and circular photogalvanic effects are studied. The 
current occurs via indirect phototransition with the participation of ripplons emission or absorp¬ 
tion. The photogalvanic tensor has strong resonances at the microwave frequency uj approaching to 
the frequencies of transitions between size-quantized subbands. The resonances are symmetric or 
antisymmetric, depending on a tensor component. Other resonances appear at cu « noJc, where n 
being integer and Wc is the cyclotron frequency. It is found that the latter resonances split to two 
peaks connected with emission or absorption of ripplons. The calculated photogalvanic coefficients 
are in accord with the experimental observed values. 

PACS numbers: 


Introduction 

The stationery surface photocurrent (in other words, 
surface photogalvanic effect (SPGE)) appears along the 
border of isotropic homogeneous bounded medium un¬ 
der the action of tilted alternative electric field [I|, [3. 
This effect was also studied in size-quantized systems [3|- 
0. Recently, SPGE attracted attention Q as a ten¬ 
tative source of the microwave-induced photoresponse 
oscillations in 2D electron gas over the liquid He sur¬ 
face (EGOHeS). The microwave-induced resistance oscil¬ 
lations (MIRO) were the subject of numerous publica¬ 
tions (see, e.g., review [ 8 Q. In particular, this effect was 
studied theoretically Q- Tl| in relation to EGOHeS (the 
theory of EGOHeS, see e.g., El). 

SPGE can be considered as another source of the ob¬ 
served photoresponse oscillations manifesting itself with¬ 
out stationary in-plane electric field. The theory of SPGE 
in EGOHeS was developed in El for the case of no mag¬ 
netic field. 

The present paper is a continuation of El with ac¬ 
counting for a strong magnetic field B directed perpen¬ 
dicular to the He surface (x, y). The forced progressive 
in-plane electron motion in a quantizing magnetic field 
is a result of the transitions between size-quantized sub¬ 
bands with synchronous directed in-plane transitions be¬ 
tween the Landau states. The mathematical reflection 
of this idea is the second order optical transition prob¬ 
ability with the participation of scattering, in particu¬ 
lar, ripplon-induced scattering. The translational mo¬ 
tion results from the interference of transition amplitudes 
caused by out- and in-plane components of alternating 
electric field. The scattering leads to a change of the in¬ 
plane electron momentum with a shift of the orbit center. 

We will mainly follow the conditions of the experiment 
0] on MIRO. We consider the electron gas of low den¬ 
sity (^ 10® cm“^) over the He3 or He4 surface. At such 
low density the electron gas is non-degenerate. The pho¬ 
ton energy is chosen close to the distance between the 


ground and the first excited size-quantized electron states 
A. The magnetic field is assumed to be weak enough so 
that the cyclotron quantum is some times less than A. 

The mentioned resonance works as a magnification fac¬ 
tor for transitions via the intermediate state. The mech¬ 
anism of SPGE can be illustrated as follows. The popu¬ 
lation of subbands changes in- or contra-phase with the 
normal component of the alternating field if the frequency 
exceeds or it is less than A/h and has ir 12 shift if fitc = A. 
This fact, together with the phase shift of the in-plane 
field component, determines the current direction. 


Problem formulation 


The phenomenology of SPGE in a magnetic field B = 
hB = B{0,0,bz) is determined by the relation for the 
current density 


j = aiRe 


(E-n(nE))(nE*) 


ra2 


aaKE - n(nE)), b](nE*) + ^[[n, [E, E*]], b], (1) 


n, [E,E* 


where n is the outer normal to the quantum well, E(t) = 
Re(Ee“*‘^‘) is the uniform microwave electric field (E = 
{Ex, Ey, Ez) is its complex amplitude). Real parameters 
ai are the functions of magnetic field value R; ai and 
as correspond to linear and 02 and 04 - to circular pho¬ 
togalvanic effects, respectively. ’’Drift” components ai 
and 02 exist in the case of zero magnetic field, while the 
’’Hall” components 03 and a\ originate from the mag¬ 
netic field action, change their signs with the magnetic 
field and vanish if B = 0. 

The current components oc 02 , <^4 can be treated sim¬ 
ilar to translational motion of a rotating wheel. Elec¬ 
tromagnetic field spin flow fc[EE*] (c is speed of light) 
transfers its angular momentum to electrons as a mo¬ 
ment of force. The friction converts this moment to the 
translational electron motion. 
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We will base on the same model of 2D EGOHeS as in 
[l^ . Electrons are attracted to He via the dielectric im¬ 
age force and the normal static electric field which com¬ 
poses Coulomb-like states xi{z) energies e;. The 
cyclotron frequency is supposed to be much lower than 
the Bohr energy. The interaction of electrons with sur¬ 
face waves (ripplons) and the homogeneous alternating 
electric field leads to the stationary surface current with 
density j. 

Eq.([T|) for current contains additive contributions oc 
and oc E^Ey. The system under consideration is 
axially-symmetric. Hence, to determine all components 
of the photogalvanic tensor, one can find only x compo¬ 
nent of the current. 

To calculate the photocurrent, we will use the approach 
first suggested by Titeika [l^. According to this ap¬ 
proach, if unpertubed electron states are localized, the 
current can be expressed via the transition probability 
between these states. Eor the x-component of the cur¬ 
rent density, one can write 




where is the probability of transitions (caused by 

perturbation) between the electron states with quantum 
numbers /3 and j3', ep and Xp are the energy and the 
center of localization, correspondingly, /(e/ 3 ) is the Fermi 
function, S is the system area, e is the electron charge. 

Let us choose the vector potential of magnetic field in 
the form of A = (0, Bx, 0) when the electron states are 
localized in the x-direction. Electron states are described 
by a set of quantum numbers /3 = (Z, n, ky), I is the num¬ 
ber of size quantized level, n is the Landau number, ky 
is the y-component of electron momentum: 

1/3) = 

VSa a 

where dimensionless oscillator functions, Xp = 

—bza^ky is the localization center (the cyclotron orbit 
center), a = \/c/\e\B is the magnetic length (we set h = 
1). These states have energies ep = £„-!-£/, where 
£„ = Wc(n -I- 1/2) is the n-th Landau level (n = 0,1, ...), 
Wc = |e|i?/TOC is the cyclotron frequency, e/ is the Z-th 
size quantization level (Z = 1, 2,.... If image attraction to 
liquid He prevails e; = — l/{2magl‘^), 

Xi(z) = 2 zexp (- 2 / 05 ) 0 ))^^^, 

X 2 iz) = z{2 - z/as) exp (-z/2aB)(2aB)“^^^, (3) 

where ub = ulme^ is the effective Bohr radius, k = 
4ki(«;i -|- K 2 )I(k 2 — Ki),Ki ^2 are dielectric constants of 
gaseous and liquid helium. Note, that the functions ([3]) 
can be used also in the presence of the normal static field, 
if to consider them as variational functions with fitting 
parameter ub- 

In our case transition probability Wp^pi is determined 
by the interaction with a microwave field and ripplons 


with the Hamiltonian TLmtit) = Uer + E{t). The Hamil¬ 
tonian of electron-ripplon interaction ZYer is 

(W)e. = 5-1/2 ^Jq(&/q-|-qq)H,(z), (4) 

q 

where 


- <lzKi{qz)), 

TTIQjq y Z 

(5) 

Jq = e*'!'’, r = {x,y), 6q, 6q are the operators of 

creation and destruction of ripplon with wave vector q 
and frequency ujq = g^/^ y^o^Jp, p is is the liquid helium 
density, (Tq is the helium surface tension coefficient. In¬ 
teraction of electron with microwave field is given by 


IP 1 

E{t) = —-b h.c. = + h.c.; (6) 

2 cj 2 

tc tc 

U = —Ev = —(E||V||-b E^x^), (7) 

UJ UJ 

where v is the operator of electron velocity. In the first 
order on the interaction Hamiltonian, the contributions 
of {U)er and E{t) to the transition probability are addi¬ 
tive and do not produce the photocurrent. Hence, the 
transition amplitude should be searched in the second 
(mixed in {lA)er and E{t)) order. In this order we have 


q,± 

^ I (^q)/ 3 / 3 'I “t“ 


“t“ CJ ^ ^g) I (■^q) 


( 8 ) 


where 



{Ici)p,PiUpi,p' 

r] -iuj + ie/ 3 i,/ 3 ' 


H/3/3i(3'q)/3i,/3' 

77 -b zw -b ze/3i./3 


iv = +0). 


(9) 


Here Iq = JqVq(z), sp^pi means ep—e'p, Nq is the ripplon 
equilibrium distribution function. Quantity {Bq)pp' is 
determined by Eq.(|9l) with change w —>■ —w and U [/+. 
For matrix elements (/q)/ 3 ,/ 3 ' one can write the following 
expressions: 


{lq)p,0' — {Jq)n,ky-,n' ,k'^{Vq)l,l', ( 10 ) 

where 


{Jq)n,ky\ri' ,ky — 

r ■\n-n'\ i{q^X+b^u,sin{2ip)/2) ib^{n'-n)if, r / 

k>X.X'-a^qyl- C e 

JnnXuq) = , (n) 

y max(n,n')! ^ minvn,rz.; 
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U^(u) is the generalized Laguerre polynomial, Uq = 
q^a'^12, kp is the polar angle of vector q = 

qcosif.qy = gsinv?); (V,);,/' = dzVq[z)xi{z)xi'{z). 
We will consider the PGE at resonance conditions, when 
the microwave frequency is close to the distance between 
size-quantization subbands I = 1 and 1 = 2. The ex¬ 
pressions for the necessary quantities QOiy, (V^) 2,2 and 
{Vq)i ,2 can be found, for example, in [l3|. 

The matrix elements of operator U are 


Ul3,l3' = —5x,X' ^1,1'+hzEy)^5n,n’ + l + 

( 'iEx bzEy')'s/Tl/Sn\n+lj ~\~ E^iVz^l.l' ( 1 ^) 


Because of low electron concentration, function fie) = 
g(A‘-e)/J’ ig the Boltzmann distribution function (/i being 
the chemical potential), / <C 1. We will consider the 
resonance case when microwave frequency is close to the 
distance between subbands with I = 1 and 1 = 2. Assum¬ 
ing that A ^ T and leaving only resonance terms one 
can reduce Ea. lfT^ to: 


^ _ nee3zi2Aa(l - e ^ 

~ r. /K 0/0 ^ 2^ ^ 


2\/2muj^iuj‘^ — uj) 


dqq^J, 


,2 j2 

n.n' 


±,n,n' 


1 , 1 


Xy/vf^iNq + IT ± -)Ei,iIm< EliujcEx - ib^ujEy) x 


Using Eas. lfTUl) . (fTTl) and (IT^ we get for PGE current: 
2e e^afb, 2'!rauJr 


Jx — 


2xaf up- 2 V 2 S 

fi^n,l)il- fien'M))Jn,n'{Nq + ^±^)qyVpi' 

q,± n,l‘,n' ,1' 


X H“ Q,l' ^ ^ ^1?) ^ 

Vpiyivz)iy,i' ivz)l,h^h,l 


Re 


f:T.{ 


rj-ioj + eiyj' ri + ioj + €iyj 


l)' 


i sig{n' — E. 


yjv! lt/n,n^ + l — 




UJc — UJ 

yjJn,n' — 1 yjTl “h 


UJc — UJ 






where E± = E^EibzEy, Sn.n' = Wc(n-n'), = ei-ei', 

sig(a:) = sgn(a: -I- 0). 

Quantities Jn,n> justify the relations: 


Jx — 


E E fien', 1 ')) 

00 

J dqq Jn,n '2 ^ ^ 

0 

T T ^q) ^ 

Vl,hivz)h,l’ iVz)l,hVlyJ' N* 


[EG 


V 77 - iw-I-eq,i' r] + iu} + ei^,i 

‘'1 

uJcEx — ibzOjEy 


7 )' 


e: 


-f (w^-a;,E^E*) U (15) 


Ul ,2 


-Tj — i5 


'—6{u}cin - n') -ojzf ojq) + 


Vi 2 

- ’—5{u}cin — n') +u}Zf uj ) + 

77 - 1-70 


- n') -STujq) 


U 2,2 


^rj — iS rj + iS 


(16) 


Here A = £2 — ei, d = A — w is the resonance detuning. 
Tie is the electron density. 

At the fulfilment of condition ojc ^ T Ea. (fT6l) can be 
simplified: 


CE 


Jx — 


ReiExE*)LOcS -I- Im(E2;E*)a;c77 - 


(13) 


6^ + rf / 

ReiEyE*)bzUJTj -|- 

where 

E = + ^)d{oJ — nu>c) + Nq9{noJc — oj) 


(17) 


/OJ \ 4/3 1 

(—) -7u”e““‘»Vip(qaB)Vi. 2 ( 903 ) 

\ oj / n\ ^ 


, (18) 


U)q — \u) — nijJc\ 


{^yjTl' H“ ^Jn.,n'-\-l y/^Jn — l,n'^^^g{j^ — y/'^qJn,n' 

i'/n/Jn,n'-l - Vn+ lJn+l,n')sigin' -n) = y/vfjJn^n' - UP) 
Using Eq. m one can rewrite Eg. (1131) in the form 


C = 


UeC^ Zl2Aaf 


12'nPa%pPiP — w^)crg^^ 


(19) 


Uq = ujg^{plao)'^^^Pl2, qas = uPj^iplao)^/^aB- 
For Vij iy) we have 

^11 (y) = 2 y^( 7 /^ - 4)“3/2 ( 7/2 - 4 ) 1/2 _ 2 arccos( 2 / 7 /) 
Ul 2 (y) = (8y2v/2/9)(47/2 _9)-5/2 [(47/2-9)1/2(9 + 87/2)- 
367/2 arccos(3/(27/)) ^ 

The comparison of Eq. o with Eq. (El) gives for the 
photogalvanic coefficients = CEap. 

UJcS UJcP 

oi = o , c-> : 02 = 


7/2 + i52 ' 
ujp 


7/2 + (52 ’ 

uj6 


“3 J^2 + ^2 ’ «4 ^2 + J2 ■ (20) 

Eg. (1751) has been obtained neglecting the Landau level 
widths. We have remained the widening 77 of the intersub- 
band distances in the prefactors only to get the finite 
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result. To include the Landau level widths, one should 
blur the delta-functions under the integral in Ea. dT^ as 
6(x) —?> In principle, widths rj and iji 

may be different, but here we will suppose rji = rj. 

Thus, the dependence of the photogalvanic coefficients 
on the magnetic field is mainly determined by the func¬ 
tion F. Fig. 1 shows this function in the case of zero elec¬ 
tron level widths. For numerical calculations we used the 
parameters of He3 and electron gas close to the conditions 
of experiment Q: = 1,K2 = 1.057, p = 0.082g/cm“^, 

(To = 0.1553erg/cm“^ and Ue = 1.4 • 10®cm“^. The pa¬ 
rameter as = 8.7 • 10“^cm is chosen to fit the intersub¬ 
band distance A = 0.38 meV corresponding to the ex¬ 
periment 0. The dependence of F{B) contains narrow 
twin peaks in the vicinity of cyclotron resonance har¬ 
monics Lo = nujc- At high temperature, the left and right 
peaks have the same amplitudes. With the drop of tem¬ 
perature, the right peak is suppressed as compared to the 
left one. The way the widening of electron levels affect 
the shape of resonances is demonstrated in Fig.2 - Fig.4. 



e, T 

FIG. 1: Function F versus magnetic field at T — 0.2K. The 
widening of the Landau levels is neglected. The double peaks 
correspond to tJ = noJc for n = 8, 7, 6, 5, 4 (from left to right). 
Insert: the fine structure of e double peak n = 5 at different 
temperatures. 


Analysis of results 

The photocurrent dependence on the parameters has 
a huge ”zoo” of resonances. First, this is the resonance 
fiw ~ A in coefficients ai (see Fa. (12(1)) ). This resonance 
originates from the participation of Ez component in 
transitions and it is not specific for the case with mag¬ 
netic field. According to our previous papers, this reso¬ 
nance appears due the intermediate state. It determines 
the symmetric or antisymmetric behavior of PGF coeffi¬ 
cients in the vicinity of resonance frequency, depending 



FIG. 2: Fig.2. The evolution of F versus magnetic field with 
the change of the Landau level width in the region of peak 
lj/loc = 5. The widening smears the twin structure. 



FIG. 3: Fig.3. The same as in Fig.2 for T = 0.03K. 

on the kind of electromagnetic field polarization (linear 
or circular). In addition to SPGF coefficients existing 
also a.t B = 0 [a, the magnetic field leads to the Hall 
photocurrent component appearance. 

Another resonance is the cyclotron resonance when the 
frequency of external field coincides with the cyclotron 
frequency. In the domain of parameters which we con¬ 
centrate here this resonance is far from our focus (and 
the experimental situation in 0)- 

One more kind of resonances is the cyclotron resonance 
harmonics w ~ nuc (corresponding to magnetic fields 
Bn = ml^c/en) on which we focus our attention (see 
Fig.l). These resonances originate from the participa¬ 
tion of the scattering in the transition processes induced 
by alternating electric field. Scattering violates the lin- 
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FIG. 4: Fig.4. The same as in Fig.2 for T = O.OIK. 

ear character of cyclotron motion resulting in cyclotron 
harmonics. Such resonances are typical for MIRO. 

According to Eq. (OH)) , the maxima (as functions of 
detuning i5) of the Hall components of the photogalvanic 
tensor a^, exceed corresponding drift components 
Qfi, 02: max5(Q;i) : max5(a2) : max5(a3) : max5(a4) = 
0.5 : 1 : n : 0.5n, where uj = nujc- This is typical situ¬ 
ation for transport in a strong magnetic field where the 
drift along the drawing force is weaker than in the Hall 
direction. The maximal values of symmetric resonances 
are 2 times larger than antisymmetric ones. 

In fact, these resonances are twinned (see insert to 
Fig.l). Note, that in quasielastic approximation the res¬ 
onances should be solitary. For example, this is the case 
when impurity scattering is take into account only. In 
the considered situation of ripplon scattering the res¬ 
onance splitting originates from the energy of emitted 
or absorbed ripplon ojq = \nujc — w|. At H —>■ i?„, 
uiq —?> 0, hence the process probability vanishes. In par¬ 
ticular, if the temperature is relatively high, T > ojq, 


F{B) oc |i? — ln(l/|H — Bn\). This explains 

deep dips between the twin peaks at nWc = to. 

On the other hand, the peaks splitting can be esti¬ 
mated as Y^CTp/ p{2n This 

estimate is consistent with insert to Fig 1. The left 
and right peaks in pairs correspond to ripplons emis¬ 
sion/absorption, accordingly. Their height ratio is {Nq -|- 
V)INq. At a T Wq the peaks heights equalize and be¬ 
have oc T ; when T —0 both peaks go down, but the left 
one tends to constant, while the right one is suppressed 
oc exp {—ojq/T). 

In the aforesaid specific conditions corresponding to 
0 the parameter C has the value C = —88 pA-cm/V^ 
at n = 5. If to suppose the level width 77 = lO^'^A 
and T = 0.2K this yields maximum of photogalvanic 
coefficient max 5 (Q; 3 ) = 2 . 5 max 5 (Q;i) = 5 max 5 (Q; 2 ) = 
2 . 5 max 5 (a 4 ) = 2.2pA-cm/V^ achieved at magnetic field 
B = 0.647T. These values are commensurable with the 
order of MIRO oscillations in Q . 

In conclusion, we have found the value of photocurrent 
along the charged liquid helium surface affected by tilted 
alternating electric field in the presence of vertical mag¬ 
netic field. Different photogalvanic coefficients represent 
responses to the linear and circular polarization of drift 
(invariant with respect to the sign of magnetic field) and 
Hall (proportional to this sign) currents. The ripplon 
scattering mechanism was taken into account. SPGE co¬ 
efficients have (symmetric or antisymmetric) their reso¬ 
nant behavior when the field frequency approaches the 
intersubband one. Besides, resonances on the cyclotron 
harmonics are observed. The current value is consistent 
with that observed in the experiment. 
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